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ABSTRACT 

We present an H -band image of the light scattered from circumstellar dust around the nearby (10 pc) 
young M star AU Microscopii (AU Mic, GJ 803, HD 197481), obtained with the Keck adaptive optics 
system. We resolve the disk both vertically and radially, tracing it over 17-60 AU from the star. Our 
AU Mic observations thus offer the possibility to probe at high spatial resolution (0.04" or 0.4 AU per 
resolution element) for morphological signatures of the debris disk on Solar-System scales. Various 
sub-structures (dust clumps and gaps) in the AU Mic disk may point to the existence of orbiting 
planets. No planets are seen in our i/-band image down to a limiting mass of 1 Mj up at >20 AU, 
although the existence of smaller planets can not be excluded from the current data. Modeling of the 
disk surface brightness distribution at i?-band and i?-band, in conjunction with the optical to sub- 
millimeter spectral energy distribution, allows us to constrain the disk geometry and the dust grain 
properties. We confirm the nearly edge-on orientation of the disk inferred from previous observations, 
and deduce an inner clearing radius <10 AU. We find evidence for a lack of small grains in the inner 
(<60 AU) disk, either as a result of primordial disk evolution, or because of destruction by Poynting- 
Robertson and/or corpuscular drag. A change in the power-law index of the surface brightness profile 
is observed near 33 AU, similar to a feature known in the profile of the f3 Pic circumstellar debris 
disk. By comparing the time scales for inter-particle collisions and Poynting-Robertson drag between 
the two systems, we argue that the breaks are linked to one of these two processes. 
Subject headings: circumstellar matter — planetary systems: protoplanetary disks — instrumentation: 
adaptive optics — stars: individual (AU Mic) — stars: low-mass, brown dwarfs 



1. INTRODUCTION 

The existence of dust disks around main-sequence stars 
has been known since the first days of t he Infrared As- 
tronom y Satellite (IRAS) mission, when lAumann et al.l 
( 1984) reported the detection of strong far-infrared 
(far-IR) excess emission toward Vega (a Lyr). Over 
200 other main-sequence stars have since been re- 
ported to possess such excesses, found almost exclusively 
with IRAS and the Infrar e d Space Observatory (e.g., 
Backnmn fc Pa rescel Il993t iMannings fc Barlowl Il998t 
Silverstond l2000t lHabi ng et al. l 120011: ISpangler et al.l 
20011; lLaureiis et alJ l2002t iDecin et alJ I2003D, though 
recently also with Spitzer (e.g.. iMever et all 120041: 
iGorlova et al.1 12004(1 . an d through ground-bas ed sub- 
millimeter observations l|Carpenter et al.l 12004^1 . Too 
old to possess remnant primordial dust, that would be 
cleared by radiation pressure and Poynting-Robertson 
(P-R) drag within several million years (Myr) in the 
absence of gas, these stars owe their far-IR excess 
to emission by "debris disks," formed by the col- 
lisional fragmentatio n of larger bodies (the so-called 
"Vega phenomenon" : iBackman fc Pa rcscc 1993, and ref- 
erences therein). Subsequent imaging at optical to 
millimeter wavelengths of the nearest sub-sample of 
Vega-like stars has resolved intricate disk-like struc- 



tures, with gaps and concent r ations (Hollan d et all 
19981 12001 iGreaves et all 119981 ISctmeider et alF " 
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Krist et al. 2000; Koerner et al. 2001; Wilncr ct al. 2002; 
Weinberger et a,Ul2002HClamnin et a,Ul2003ft . 

The most favored explanation for such structures 
is the gravitational perturbation by embedded plan- 
ets orbiting at semi-major axes comparable to the disk 
size llMoro-Martm fc Malhotral2002l:lKenvon fc Bromlevl 
2004). The existence of perturbing planets may be re- 
vealed by clumps of dust trapped in mea n motion res- 
onance s, as ha s been suggested f or Vega (iWilner et all 

20021 IWvattJ |20?)1 e Eri JOzernov et all 1200(1 

Quillen fc Thorndikd 120021) . and Foma lhaut (a PsA; 
Wvatt fc Dentll2002l: iHolland et al.ll2003ft . and observed 
by the Cosmic B ackground Explor er satellite along 
the Earth's orbit (Re ach et alJ 11995^) . Stochastic col- 
lisions between large planetesimals that result in dust 
clumps lasting several hundreds of orbital perio ds are 
anoth er means of producing disk asymmetries (jSternl 
1996). Spiral density waves (as seen in the disk of 
HD 141569A, and inferre d around /3 Pic: iClamoin et all 
20031 iKalas et all l2000(l . and warp s in the disk in- 
clination (as in the disk of (3 Pic; IHeap et all 120001 
IWahhai et all l200l. may indicate perturbation by 
nearby stars l|Ka las el; alpQQOLj^Ql UKenvon fc Bromlevl 
120021 lAugereau fc Papaloizoull2004|) . Finally, dust mi- 
gration in a gas-rich disk can produce azimuthally sym- 
metric structures, as observed in the HR 4796A circum- 
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stellar disk ijTakeuchi fc Artvmowiczll2001[) . 

High-resol ution imaging observat ions, suc h as those of 
HR 4 796A Schneider et all Il99l. Pic fHeap et al.l 
12000^ . HP 141569A llWeinberger et alJ 119991) and 
TW Hya (|Krist et all200ft TWcinbcr ger et al]l2002t) with 
the Hubble Space Telescope (HST) can help single out 
the most likely physical process behind the disk mor- 
phology. The resolution achievable with adaptive optics 
in the near-infrared on large ground-based telescopes ri- 
vals that of HST, and is the method employed in this 
paper for investigating disk structure. 

The young (8-20 Myr; iBarrado v N ayascu^s et al.l 
Il999t IZuckerman et al]l200ll) Ml V l|Keenanlll983l) star 
AU Mic has a know n 60/zm excess fr om IRAS, likely 
due to orbiting dust l)Song et al.ll2002|) . Because of its 
relative proximity (Hi pparcos distance of 9.94±0.13 pc; 
iPerrvman et alJ Il997|) . AU Mic is a good target for 
high-resolution imaging of scattered light to characterize 
the circumstellar disk morphology . Rec ent 450/im and 
850/im observations by iLiu et all ((21)04) confirmed the 
existence of dust, and follow-up i?-band (0.65/im) coron- 
agraphic imag ing revealed a nearly edge-on disk extend- 
ing 210 AU (|Kalas et alJl2004j) from the star. Because 
the age of AU Mic is larger than the collision timescale 
between pa r ticles in the disk (0.5-5 Myr at 200 AU), 
iKalas et all 1)20041) infer that most of the dust particles 
have undergone at least one (destructive) collision, and 
hence the AU Mic disk is a debris disk. However, because 
the P-R time scale for 0.1-10^m particles at >100 AU 
from the star is greater than the stellar age. IKalas et alJ 
expect that most of the disk at large radii consists of 
primordial material. ILiu et all find a fractional infrared 
luminosity, L[r/L* = 6 x 10~ 4 and fit the far-IR to sub- 
millimeter excess by a 40 K modified blackbody with 
constant emissivity for A < 100/xm and following A~ - 8 
for longer wavelengths. From the lack of excess at 25/im, 
ILiu et all infer an inner disk edge at 17 AU from the star, 
or 1.7" at the distance of AU Mic. They speculate that 
such a gap may have been opened by an orbiting planet, 
that, given the youth of the system, could be detectable 
in deep adaptive optics (AO) observations in the near 
IR. Such high contrast observations could also be used 
to search for signatures of planet/disk interaction. 

The AU Mic circumstellar disk is not resolved with 
the 14" b eam of the JCMT/ SCUBA observations of 
ILiu et ail , and the IKalas et al.1 optical coronagraphic ob- 
servations are insensitive to the disk at separations <5" 
because of the large sizes of their occulting spots (di- 
ameters of 6.5" and 9.5"), and because of point-spread 
function (PSF) artifacts. Taking advantage of the higher 
angular resolution and dynamic ra nge achiev able with 
adaptive optics on large telescopes, iLiul l)2004[) used the 
Keck AO system to investigate the disk morphology at 
separations as small as 15-20 AU from the star. We 
pres ent our own set of Keck AO data that confirms 
iLiuf s observations, and places upper limits on the pres- 
ence of potential planetary companions. In addition, we 
combine our spatially resolved ff-band information with 
the i?-band imaging data from IKalas et al.l l)2004j) . and 
with the optical to sub-millimeter spectral e nergy dis- 
tribution (SED) of AU Mic from ILiu et alJ l|2004j) . to 
put self-consistent constraints on the disk morphology 
and the dust properties, as done prev iously for Pic 
ijArtvmowicz. Burrows, fc Parescelll989D . We use a full 



three-dimensional radiative transfer code to model simul- 
taneously the SED and the ff-band and i?-band surface 
brightness profiles (SBPs). We find that our model con- 
sisting of a single dust population does not reproduce 
the observed break in the 7J-band SBP, whereas a two- 
component dust model, as proposed for Pic, fits the 
data well. Drawing from a comparison with the Pic 
system, we deduce that dynamical evolution of the disk 
provides the simplest explanation for the morphology of 
the SBPs of both disks. 

2. OBSERVATIONS AND DATA REDUCTION 

We observed AU Mic at H band (1.63/im) with the 
NIRC2 inst rument (Matthews et a l., in prep.) and the 
AO system l)Wizinowich et al1l2000j) on the Keck II tele- 
scope. The data were acquired on 5 June, 2004 un- 
der photometric conditions. We employed coronagraphic 
spots of different sizes (0.6"-2.0" diameter) to block out 
the light from the star. The observations were taken with 
the wide, 0.04" pix -1 camera in NIRC2, which delivers 
a 41" x 41" field of view (FOV) on the 1024x 1024 InSb 
Alladin-3 array. 

We obtained nine 54 sec exposures at H band, three 
with each of the 0.6", 1.0", and 2.0"-diameter corona- 
graphic spots. We observed the nearby (2.4° separation) 
M2/3 III star HD 195720 as a PSF standard, with simi- 
lar colors, but 0.9 mag brighter than AU Mic at H . We 
spent equal amounts of time on the target, on the PSF 
star, and on sky. The observations were carried out ac- 
cording to the following sequence, repeated 3 times: 3 
exposures of AU Mic, 3 exposures of HD 195720, and 3 
exposures of the blank sky (taken at three dithered posi- 
tions: 60", 50", and 40" away from the PSF star). The 
total on-source exposure time was 8.1 min. Throughout 
the observations, the field rotator was set in "position 
angle mode," preserving the orientation of the sky on 
the detector. The image quality was estimated from the 
Strehl ratios of point sources observed at higher spatial 
resolution (with the narrow camera, 0.01"/pix) at the 
beginning of each night. Our H band images had Strehls 
of 17-20%. 

Data reduction followed the standard steps of sky- 
subtraction, flat-fielding, and bad-pixel correction. The 
images in the individual sets of 3 exposures were then 
median-combined to improve the sensitivity to faint ob- 
jects. The AU Mic dust disk was barely discernible at 
this stage (Figure^,). To enhance the visibility of the 
dust disk, we subtracted the stellar PSF. The PSF was 
obtained by first rotating the image of HD 195720 to 
match the orientation of the diffraction pattern in the 
AU Mic image, and then scaling it by a centrally sym- 
metric function f(r) to match the radial dependence of 
the AU Mic profile. The objective of the scaling was 
to compensate for variations in the seeing halo between 
AU Mic and the control star caused by changing atmo- 
spheric conditions and fluctuating quality of the AO cor- 
rection. The function f(r) was obtained as the median 
radial profile of the ratio of the AU Mic to the HD 195720 
images, with the telescope spikes and the edge-on disk 
masked. Figure [21 shows f(r) for the images with the 
three different coronagraphic spot sizes. Two remarks 
on this procedure should be made here. First, the func- 
tion f(r) does not vary by more than 15% for any given 
spot size, and for large radii tends to 0.4 — the -ff-band 
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flux ratio between AU Mic and the PSF star. Second, 
because the AU Mic disk happens to be edge-on, the cen- 
trally symmetric scaling of the PSF does not introduce 
any spurious features in the result, and so does not in- 
terfere with the morphology of the disk. The procedure 
would not be viable for disks that are far from edge-on. 

The scaled version of HD 195720 for the image with 
the 0.6" coronagraph is presented in Figure Panel 
(c) of Figure n shows the final AU Mic image, obtained 
by median-combining all coronagraphic exposures. In 
panel (d) a digital mask has been employed to enhance 
the appearance of the circumstellar disk. The mask en- 
compasses the innermost 1.7" from AU Mic, as well as 
the hexagonal diffraction spikes from 2 of the 3 image 
sequences where they did not subtract well (as seen in 
panel (c) of the Figure). 

For flux calibration we adopted the 2MASS magni- 
tude for AU Mic (H = 4.831 ± 0.016) and relied on the 
residual transmission of the NIRC2 coronagraphic spots. 
The PSF star, HD 195720, is saturated in 2MASS, and 
therefore unusable for flux calibration. We measured the 
flux from AU Mic through the 1.0" and 2.0" corona- 
graphs in a 6 pix (0.24") diameter aperture. From non- 
coronagraphic images taken with the 10 mas/pix cam- 
era on NIRC2, we found that this aperture contained 
68% of the total power in the PSF. The -ff-band trans- 
missivity of the 2" NIRC2 coronagraph was measured 
at (7.32 ± 0.24) x 10" 4 (extinction of 7.84 ± 0.03 mag). 
AU Mic is known to exhibit a large U-band photometric 
amplitud e (0.35 mag) due to star spot s with a period of 
4.9 days l)Torres fc Ferraz Mellolll973|) . Although we do 
not expect measurable variability over the «1 hour time 
span of our observations, our absolute flux calibration is 
uncertain. Nevertheless, in the near IR the contrast be- 
tween the spots and the stellar photosphere is less pro- 
nounced than in the optical, so the uncertainty is also 
smaller: likely of order <0.1 mag. 

The absolute orientation of the dust disk arms was cal- 
ibrated through observations of a binary star standard, 
WDS 18 055+0230, with well-known orbital ele ments 
(grade 1: lHartkopf & Masonll2003t iPourbaixl EoOOh The 
y-axis of the NIRC2 detector was measured to be off- 
set by 1.24° ± 0.10° clockwise from North. All position 
angles quoted below have been corrected for this offset. 

3. RESULTS & ANALYSIS 

3.1. Circumstellar Dust Morphology 

The disk is seen out to a distance of ~6" (60 AU) from 
the star in our combined ii-band image (Figure ^^d). 
Inwards it can be traced inwards to « 1.7" (17 AU) from 
the star, at which point residual speckle noise from the 
PSF subtraction overwhelms the emission from the disk. 
Thus, our imaging data cannot directly test the exis- 
tence of the p roposed disk cle aring interior to 17 AU 
from the star ijLiu et al.ll2004j) . We confirm the sharp 
mid-p lane morphology of the disk ijKalas et al.ll2004l iLiul 
2004), indicating a nearly edge-on orientation, and re- 
solve the disk thickness, with the SE arm appearing 
somewhat thicker (FWHM = 2.8-4.4 AU) than the NW 
arm (FWHM = 2.2-4.0 AU). There is also evidence of an 
increase in the FWHM of each of the arms with separa- 
tion: from 2.2-2.8 AU at 20 AU to 4.0-4.4 AU at 40 AU 
from the star, indicating a potential non-zero opening 
angle of the disk. Within 5" of AU Mic, the position 



angles (PAs) of the two sides of the disk are nearly 180° 
away from each other: we measure PA = 310.1°±0.2° 
for the NW arm, and PA =129.5°±0.4° for the SE arm. 
These PAs are in agreement with those reported in iLhJ 
l|2004L 311.4°±1.0° and 129.3°±0.8°), though more ac- 
curate, likely as a result of our proper calibration of the 
orientation of the NIRC2 detector (Section 2). 

The radial SBP of the disk was measured on the 
reduced image (Figure before applying the digital 
mask) using the IRAF task polyphot. The photom- 
etry regions are indicated by the rectangles overlapped 
onto a contour map of the image in Figure Et- We used 
4 pix x 12 pix (0.16" x 0.48") rectangular apertures, 
where the long side of the rectangular regions was cho- 
sen to span 1-2 FWHMs of the disk thickness, and was 
aligned normally to the disk arm. The distance between 
the aperture centers was 4 pix (0.16"). Even though the 
image was PSF- and sky-subtracted, to offset for flux 
biases introduced by the centrally symmetric scaling of 
the PSF (Section 2), we employed additional background 
subtraction, with the background flux estimated as the 
median pixel value in 0.1 6"- wide concentric annuli cen- 
tered on the star. The photometric uncertainty was esti- 
mated as the quadrature sum of the standard deviation 
of the background and the photon noise from the disk. 
For the standard deviation of the background we adopted 
the root mean square (r.m.s.) of the pixel values in the 
annulus, multiplied by \Ztt/2 (to pr operly account for 
the s tandard deviation of the median; .Kend all fc Stuart! 
Il977|) . and normalized by the size of the photometry 
aperture. 

The radial SBPs of the NW (upward pointing trian- 
gles) and SE (downward pointing triangles) arms of the 
projected disk are shown in Figure 0] The two SBPs 
agree well throughout the region over which we can trace 
the disk (17-60 AU). Unlike the observed shape of the R- 
band SBP at 50-210 AU, the 7J-band SBP at 17-60 AU 
from the star cannot be fit by a single power law. Instead, 
the SBPs of both the NW and the SE arms "kink" and 
flatten inwards of 30-40 AU, with the transition being 
more abrupt in the SE arm at «33 AU, and more grad- 
ual in the NW arm. While the power-law exponent of the 
mean SBP over the entire range (17-60 AU) is — 2.3±0.2, 
over 17-33 AU separations it flattens to — 1.2±0.3, while 
over 33-60 AU, it increases to —4.0 ±0.6. These are con- 
sistent with the measurements of ILiul ()2004f) over the 
same separation range. 

A closer look at the SBP of the NW and SE arms re- 
veals several small-scale asymmetries, all of which can 
be linked to regions of non-uniform brightness in the 
AU Mic disk (Figure |3]i). The sub-structure is enhanced 
by scaling the reduced image by a radially symmetric 
function centered on the star, with magnitude propor- 
tional to the radius squared (Figure |3b) . The lettered 
struct ures denote features identified in the (deeper) im- 
age of lLiul (|2004l) : clumps of enhanced emission (A, and 
C), a gap (B), and a region elevated with respect to the 
inner disk mid-plane (D). In general, we confirm the pres- 
ence of these features in the AU Mic disk, although the 
gap (B) and the clump (A) appear misplaced by ~5 AU 
toward the star in our image. The NW arm also looks 
more uniform in brightnes s between 17-40 AU in our 
image compared to that in lLiiJ l)2004l Figures 3 and 4), 
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where clump A is very prominent. These discrepancies 
may be caused by residual speckle noise from the PSF 
subtraction i n eit her data set. In addition to the features 
described bv iLiul . we see a faint clump in the NW arm 
at >50 AU from the star, coincident with the location of 
the bump in the SBP of th is ar m (Figure^}. The concen- 
tration is not reported bv ILiul and being at a relatively 
low signal-to-noise (~3), may be a noise spike. 

3.2. Disk Luminosity, Optical Depth and Geometry 

The integrated disk brightness (over 17-60 AU from 
the star) is 2.7±0.8 mJy at H band; hence L sca t/L* = 
2.3 x 10~ 4 . This is comparable to the fractional dust 
luminosity in emitte d mid-IR t o sub - mm light, fd = 
Lir/L* = 6 x 10~ 4 iLiu et all l)2004|) . and hence sug- 
gests that the disk mid-plane may be optically thin to 
radiation at wavelengths as short as ~1.5/im, at the 
peak of the AU Mic spectrum. Indeed, fd is simi- 
lar to that of other known debris disks (10~ 5 — 10~ 3 ; 
e.g.. iSvlvester fc Mannings! 120001 lHabing et all 120011; 
Soan gler et al.ll200ll) . all of which are optically thin to 
ultra-violet and optical light in the direction perpendic- 
ular to the disk plane (t_l < 1). In the mid- plane, the 
optical depth of grains along a radial line from the star to 
infinity is tii ~ fd/ sin<5 if the grains are in a "wedge" or 
"flaring" disk wit h thickness proportional t o radius and 
opening angle 28 l)Backman fc Parescell993lh Because of 
the generally unknown viewing geometry of circumstel- 
lar disks, ry tends to be poorly constrained. Assuming 
edge-on orientation (i = 90°), we can estimate the max- 
imum allowed opening angle 26 from the observed disk 
thickness. For smaller values of i, 5 will be smaller be- 
cause of projection effects. Assuming a perfectly flat, 
thin disk, we find a lower limit on the inclination i > 87° 
over 20-50 AU. The projected appearance of an inclined 
disk of zero thickness would however be inconsistent with 
the apparent thickening of the disk with increasing sep- 
aration fSection l3.ll) . Therefore, the disk likely has non- 
zero scale height and/or opening angle, and is viewed 
within only a degree of edge-on. To put a limit on 5, 
we observe that at a radius of 40 AU the disk FWHM is 
~4 AU. Thus we obtain that for i < 90°, S < 3°. Hence, 
tii > 6 x 10~ 3 , and probably less than unity, i.e., the disk 
is optically thin in the radial direction. 

3.3. Detection Limits on Sub-Stellar Companions 

Dynamical influence by embedded planets is a fre- 
quently invoked explanation for substructure in dust 
disks. Because of its youth, proximity, and late spec- 
tral type, AU Mic is an ideal target for direct imaging of 
planets. However, no point sources are seen in our com- 
bined 8.1 min PSF-subtracted H band exposure. Fig- 
ure delineates our 5a sensitivity limits as a function of 
angular separation from the star in the PSF-subtracted 
image. The la level at each distance is defined as the 
r.m.s. scatter of the pixel values in one-pixel wide an- 
nuli centered on the star. This was divided by V28.3 — 1 
to adjust for the finite size of the aperture used for pho- 
tometry (six-pixel diameter, or an area of 28.3 pix 2 ). The 
sensitivity to point sources within the disk is up to 1 mag 
poorer because of the higher photon noise. This is shown 
in FigureJSJwith the points, each of which has the photon 
noise from the disk signal added in quadrature. These 



detection limits, calculated in a statistical manner, were 
confirmed through limited experiments with artificially 
planted stars. 

At the location of the inferred gap in the SE arm (~25- 
30 AU), we can detect planets down to 1 Jupiter mass 
(Mj up ) for an assumed stellar ag e of 10 Myr and usin g 
brown-dwarf cooling models from iBurrows 
Dynamical models of planet-disk interactions in other 
systems exhibiting similar di sk morphology (e Eri, a Lyr) 
require planets 0.1 -3 M Jup ijQuillen fc Thorndikd 12002: 
iWilner et alJ 12002) . Provided that the clumps in the 
AU Mic disk are caused by such a planet, our point source 
detection limits constrain its mass to the lower part of 
this range. 

A faint candidate companion is seen around our PSF 
star, HD 195720 (H = 3.88 ± 0.24 from 2MASS). The 
object is 9.5±0.2 mag fainter at H, at a projected sep- 
aration of 1.19" and PA of 81° (Figured). Given that 
HD 195720 is a distant giant star (spectral type M2-3 III 
from SIMBAD), if associated, the companion would be 
a main-sequence K star. Because of the large magni- 
tude difference, the presence of the projected companion 
does not affect our PSF subtraction or the analysis of the 
AU Mic circumstellar disk. 

4. DUST DISK MODELING 

It has already been suggested that the structures in 
the AU Mic disk (dust concentrations, gaps, vertically 
displaced clumps) are likely signposts of the existence 
of perturb ing planetary-mass bodies in the AU Mic disk 
l)Liull200"H). The proposed clearing in the disk inwards 
of 17 AU l)Liu et alJl2004|) supports such a hypothesis. 
From our imaging data we cannot trace the disk to sep- 
arations <17 AU to directly test the existence of a gap. 
However, by combining our knowledge of the optical to 
sub-mm data on the AU Mic debris disk with an ap- 
propriate model, we can still probe some of the physi- 
cal properties of the disk, including the size of the inner 
gap. In this section we present results from a three- 
dimensional continuum radiative tr ansfer code, MC3D 
(Wolf fc HenningJ 120001 iWolfl l2003|) . to simultaneously 
model the AU Mic SED and the scattered light from 
the disk, and to place constraints on the dust grain size 
distribution, the radial particle density distribution, and 
the inner disk radius. 

4.1. Model and Method 

The MC3D code is based on the Monte Carlo 
method and solves the radiative transfer problem self- 
consistently. It estimates the spatial temperature distri- 
bution of circumstellar dust, and takes into account ab- 
sorption and multiple scattering events. Given the non- 
vanishing mid-plane optical depth of the AU Mic disk 
fSection l3.2fl . we believe that the use of a multi-scattering 
approach is warrant ed. The code employs the concept of 
enforced scattering ijCashwell fc Everettlll959j) . where in 
a medium of optical depth r, a fraction e~ r of each pho- 
ton leaves the model space without interaction, while 
the remaining part (1 — e~ r ) is scattered. The code 
is therefore applicable to the low-density environments 
of circumstellar debris disks. The dust grains are as- 
sumed to be spher ical with a power-law size distribu - 
tion, n(a) oc a~ 3,5 (Mathis. Rumpl. fc Nordsieckl ll977). 
We used a standard inter-stellar medium (ISM) mixture 
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of 62.5% astronomical silicate and 25% + 12.5% graphit e 
Ipraine fc Malhot.r Jffl53 IWeingartner fc prainflkOOl I 
with optical properties from iDraine &; Lea 0984). The 
extinction and scattering cross sections, and the scatter- 
ing distribution functio n are modeled following t he Mie 
scattering algorithm of Bohrcn & Huffman (1983). 

We use the MC3D code to model the R and H 
band scattered light in the AU Mic disk, and the IR 
to sub- mm excess in the SED. The AU Mic photo- 
sphere is best approximated by a 36 00 K NextGen model 
(Hauschildt, Allard Tfc Baronll999|) . as expected from its 
spectral type (Ml V). The fit was performed over the 
1-12/im wavelength range, where the emission is photo- 
spheric. Data from the literature at shorter wavelengths 
were ignored, as they are not taken simultaneously, and 
hence are strongly affected by the V = 0.35 mag vari- 
ability of the star. By matching the model i^g-band flux 
density to that of a blackbody of the same temperature, 
and adopting the Hipparcos distance of 9.94 pc to the 
star, we find that its luminosity and radius are 0.13L Q 
and O.93i?0, respectively. For the debris disk we adopt 
a flat (unflared) geometry with a number density pro- 
file proportional to r~ 7 , where r denotes radial distance 
from the star, and 7 is a constant. We set the outer ra- 
dius of the model to 1000 AU, so that it is larger than the 
size of the i?-band scattered light emission (210 AU), and 
of the JCMT/SCUBA beam used for the sub- mm mea- 
surements (FWHM of 14" = 140 AU). The disk inclina- 
tion and opening angle were already constrained in Sec- 
tion For our modeling purposes we assume i = 89°, 
(5 = 0°, and a flat disk model with a constant scale height 
h=0.8 AU. We find that models based on these parame- 
ters approximate the mean observed disk thickness well. 

The remaining free parameters in the disk model are 
the exponent of the volume density profile 7, the dust 
mass Mdust, the minimum and maximum dust grain sizes 
a m in and a max , and the inner radius r in . A fit to the 
mean NW and SE SBP between 17-60 AU results in a 
best-fit power-law index of — v = —2.3 ±0.2, indicat- 
ing that the number density profile varies as r - 1 - 3±0 - 2 j 
i.e., 7 = v — 1 = 1.3 ± 0.2 (as is tru e for an edge- 
on disk of isotropically scattering gra ins ; INakanol 1*1990: 
IBackman. Witteborn. fc Gillettjll992|) 1 . This value is in 
agreement with the rang e inferred for P-R drag domi- 
nated disks (1.0-1.3; e.g.. lBriggsllT96^ . G iven the error 
on our fit, we decide to fix the value of the power-law 
index at the theoretically expected value of 7 = 1.0 
for a continuously re plenished dust cloud in equilib- 
rium under P-R drag llLeinert. Roser. fc BuitragolH983t 
Bac kman fc GilletHll987ir The effects of varying 7 are 
considered at the end of Section 14.21 

4.2. Breaking Degeneracies in the Model Parameters 

We subsequently follow a trial-and-error by-eye opti- 
mization scheme to determine the values of Md us t, amin, 
o mM , and r- m . With a sophisticated dust disk model 
containing many parameter choices, it is possible to find 
combinations of parameters that have degenerate effects 
on the SED and/or on the SBP. By fitting simultaneously 

Forward scattering, to the extend to which it is characteristic 
of the dust grains in the AU Mic disk, tends to increase 7. Even 
though forward scattering is ignored in the approximation 7 = v — 
1, it is modeled by the MC3D code, where its amount is determined 
by the input grain parameters and Mie theory. 



the SED and the imaging data we can avoid some, but 
not all, of the complications. Here we discuss the specific 
degeneracies and how we can break them via the observa- 
tional constraints in hand. We first consider the interac- 
tion between Md us t, Omin, and a max , which are strongly 
degenerate. We then consider the effect of changing r- m , 
which is more weakly coupled with the rest of the pa- 
rameters. Finally, we extend our discussion to consider 
variations in the power-law index 7, which is otherwise 
kept fixed during the modeling. 

The dust mass Mdust, and the minimum and maxi- 
mum grain sizes, a mm and a max , have degenerate effects 
on both the SBP and the SED. Decreasing M^ ust , or in- 
creasing Omin or a max , result in a decrease in the amount 
of mass residing in small grains (the bulk of the scat- 
terers), and lowers the flux of the SBP. Each of these 
changes similarly lowers the thermal SED flux. However, 
dust mass variations can be disentangled from grain size 
variations because of the different magnitudes of their ef- 
fects on the SBP and on the SED. Optically thin thermal 
emission is a more accurate proxy of dust mass, whereas 
optical-IR scattering is more sensitive to small differences 
in the mean grain size. We therefore constrain Mdust 
from the sub-mm data, while we use the color and ab- 
solute flux of the scattered light to determine a m i n and 
Umax- Here we should note that Md us t represents only 
the dust mass contained in grains comparable in size or 
smaller than the maximum wavelength (A max ) at which 
thermal emission is observed. In the case of AU Mic, the 
current ly existing longe st-wavelength observations are at 
850/im (|Liu et all2004D . Consequently, we are free to ad- 
just a min and a max , as long as a min < a max < 1 mm. We 
will not consider cases for which a max > 1 mm. 

As a first step in finding the optim um model parame- 
ters, we confirm the lLiu et alJ l)2004[) estimate of the dust 
mass, Mdust = 0.011M®, calculated from the 850/im flux. 
This value matches the sub-mm data points for a wide 
range (an order of magnitude) of minimum and max- 
imum grain sizes, whereas changing Md us t by a factor 
of >1.5 introduces significant discrepancies from the ob- 
served 850/im emission. As a next step, we constrain 
Omin by modeling the optical-near-IR color of the scat- 
tered light in the overlap region (50-60 AU) between the 
iKalas et alJ i?-band and our H-b&nd data. Because the 
size of the smallest grains is likely comparable to the 
central wavelengths of the R and H bands, R — H is 
a sensitive diagnostic for a nnn . We smooth the R- and 
-f/-band model images to the respective image resolutions 
(1.1" at R and 0.04" at H), and use the ap propriate aper- 
ture widths (1.2" at R |Kalas et al.ll2004] and 0.48" [this 
work]). For the adopted grain size distribution we find 
that models with a m ; n « 0.5/im, with a probable range 
of 0.3-1.0/im, best approximate the disk color. Having 
constrained Md us t and a m ; n , we find that a max = 300/im 
matches best the R- and H-b&nd flux levels of the disk. 
The probable range for a max is 100-1000/im. Since we 
have not considered models with a max > 1 mm, we can- 
not put an upper limit to the maximum grain size. 

The inner radius rj n of the disk is degenerate with the 
mean dust grain size in the SED. Greater values of r; n 
decrease the mid-IR flux and shift the peak of the excess 
to longer wavelengths (Figure Et), as do greater values 
of a m i n (Figure EJd) and a max . Because the AU Mic disk 
is optically thin in the mid-plane fSection l3.2|) . the inner 
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disk radius has no effect on the flux and color of the SBP, 
facilitating the isolation of the n n parameter. Having 
already determined a m in and a max , we find r- ln s» 10 AU: 
smaller than the 17 AU gap estimated f rom t he single- 
temperature blackbody fit in iLiu et alJ (|2004f) . A firm 
upper limit of r- m < 17 AU can be set based on the fact 
that we do not observe a decrease in the intensity of the 
scattered light with decreasing separatio n down to 17 AU 
(Figure El see also Figure 2 in lLiull2004|) . 

Finally, although we had fixed the value of the power- 
law index 7, a brief discussion of its variation is war- 
ranted given the change in the SBP with radius. The 
action of 7 on the SED is degenerate with the dust size, 
and with the radius of the inner gap. Larger values of 7 
are degenerate with smaller particles and smaller inner 
gap radii (Figure EJ). Given that at <33 AU the power- 
law index of the SBP decreases to —1.2 (i.e., 7 w 0.2), 
and that most of the mid-IR flux is produced close (10- 
20 AU) to AU Mic, the inner disk clearing may there- 
fore be smaller than 10 AU in radius. Indeed , recent 
HST scattered light imaging (|Krist et al J 12004(1 detects 
the disk in to 7.5 AU (although the authors invoke for- 
ward scattering to account for the apparent filling in of 
the inner gap). Compounding this with evidence for an 
increasing minimum grain size with decreasing radial sep- 
aration fSection l5.1jl . we find that inner gap sizes as small 
as r in ~ 1 AU cannot be ruled out. 

Table ^ lists our preferred model parameters for the 
AU Mic star-disk system. The optical depth of the model 
along the disk mid-plane is tii » 0.08 at both R and H 
bands, in agreement with the estimates in Section 13.21 
The model SBP and SED are over-plotted on the data 
in Figures 01 and respectively. Figure 0t shows the 
noiseless scattered light model of the AU Mic disk with 
the same greyscale and contour spacing as the image in 
Figure Eli. The model disk extends to larger radial sep- 
arations than the AU Mic disk: an effect of the steeper 
power law of the AU Mic SBP at >33 AU. 

5. DISCUSSION 

Two important new results are evident from our simul- 
taneous modeling of the SBP and SED of the AU Mic 
debris disk: (1) there is a pronounced lack of small 
(<0.3^m) grains in the inner disk, and (2) the radius 
of the inner clearing may be smaller (1-10 AU) than es- 
timated JT7_AU) from a simple black-body fit to the IR 
excess (|Liu et al.l [2 004). The latter point was already 
discussed in Section l4~2l and, given the shortcomings of 
our model in reproducing the changing slope of the SBP, 
will not be belabored further. Here we discuss the de- 
rived minimum grain size along with recent evidence for 
its dependence on disk radius. We then focus on the 
change in slope of the SBP of the AU Mic debris disk, 
and draw a parallel with the ft Pic system. We propose 
that identical dynamical processes in the two debris disks 
can explain the observed homology. 

5.1. Minimum Grain Size as a Function of Disk Radius 

From the model fits to the color and absolute flux of 
the scattered light from the AU Mic debris disk, we find 
that the dust grains are between a m i n = 0.5^Q 2/ini and 

Qmax = 300^2ooM m i n s i ze (although the 1 mm upper 
limit on the maximum grain size is not robust). In reality, 



our constraints on the grain parameters are valid only 
over the 50-60 AU region, where we have information 
from both the R — H color of the scattered light and the 
SED. We have very few constraints for the outer disk 
(> 60 AU), which is seen only in i?-band scattered light 
and is too cold to be detected in emission at wavelengths 
<1 mm. 

Shortly before receiving the referee report for this 
paper, sensitive high-resolution 0.4 -0. 8 um images o f 
AU Mic became available from HST (|Krist et alJl2004(l . 
These show the debris disk over 7.5-150 AU separations 
from the star and thus provide complete overlap with our 
AO data. A brief discussion of the two data sets in the 
context of the minimum grain size is therefore warranted. 
For consistency with iKrist et all we re-did our surface 
photometry with the 0.25" x 0.25" apertures used by 
these authors. Because the two data sets have similar an- 
gular resolution (0.04" vs. 0.07"), and because the aper- 
ture size is much larger than the FWHM of the PSFs, the 
systematics of the photometry should be negligible. The 
IKrist et al] HST F6 06W (0.59/xm) data are consistent 
with the lKalas et alJ ground-based i?-band (0.65/im) ob- 
servations over 50-60 AU from the star, and thus our con- 
clusions about the grain sizes in this region remain un- 
changed. However, a comparison of the iJ-band and the 
i 7 '606IU-band SBP over the region of overlap shows that 
the F606W^ — H color changes from 2.9±0.2 mag (i.e., 
approximately neutral, since R — H = 2.9 for AU Mic) 
at 17-20 AU to 2.0±0.3 mag at 50-60 AU. That is, the 
debris disk becomes increasingly bluer at larger radii . 
The effect is gradual and is also reported in IKrist et all 
where it is observed over a narrower wavelength range 
(0.4-0.8/xm) at 30-60 AU from the star. 

The neutral color of the dust at 20 AU indicates that 
the majority of scatterers there are larger than 1.6/im. 
Compared with the minimum grain size (0.5^2^^) that 
we derived at 50-60 AU, this indicates that at smaller 
separations grains are bigger. Such a dependence of 
grain size on radius would further imply that the ra- 
dius at which the SBP changes power-law index should 
be wavelength-dependent, occurring farther away from 
the star at shorter wavelengths. Evidence for this may 
indeed be inferred from a comparison between the HST 
and Keck AO data: in the F606W ACS filter the break in 
the SBP is seen at «43 AU (|Krist et all2004H . whereas at 
iJ-band it occurs near 33-35 AU ( Section EHl lLhjff2 0041 . 
The indication that the minimum grain size decreases 
with disk radius is thus confirmed from two independent 
observations. 

Particles smaller than a m i n may be removed as a result 
of either coagulation into larger particles (grain growth) , 
destruction by P-R and/or corpuscular drag, or radiation 
pressure blow-out. Given that a m ; n is larger than the 
radiation pressure blow-out size (0.14/zm, for a radiation 
pressure to gravity ratio /3=0.5 and a grain density of 
2.5 g cm -3 ), grain collisions and drag forces dominate the 
dynamics of >0.14/im grains around AU Mic. Therefore, 
the origin of the sub-micron grains scattering visible light 
at wide separations (where the collision and P-R time 
scales are longer than the age of the star) may also be 
primordial: rather than being blown out from the inner 
disk, these grains may be remnant from the proto-stellar 
cloud that never coagulated beyond an ISM grain size 
distribution. 



AO Imaging of AU Mic 



7 



5.2. The Change in the SBP Power-Law Index: a 
Comparison with Pic 

It is not surprising that our preferred model cannot 
reproduce the detailed structure of our high angular res- 
olution IR image. The model parameters were found 
only after a coarse sampling of the parameter space, and 
through a number of simplistic assumptions that merely 
approximate the physical conditions in the AU Mic debris 
disk. In particular, under the assumption of a uniform 
grain size and density distribution over 10-1000 AU, the 
MC3D model cannot mimic the SBP slope change at 
~33 AU an d the clump y substructure over 17-60 AU 
described in iLiul l)2004|) and confirmed in Section 13.11 
While the dust clumps are high-order perturbations that 
may require dynamical considerations for proper model- 
ing, the change in the SBP potentially could be explained 
in the framework of existing dust disk scenarios. 

The occurrence of the power-law break at similar radii 
in the SBPs both arms of the projected disk suggests 
that this is a ring-like structure surrounding the star, 
rather than a discrete feature at one location in the disk. 
Models involving dynamical interaction with planets 
have been proposed to explain ring-like structures in cir- 
cumst ellar disks (e.g . . IRoques et al.lll994t iLiou fc Zoqkl 
Il999t iKenvon et all 119991: iKenvon fc BromlevT boO^. 
and the clumpy structure of the AU Mi c disk doe s 
suggest the presence of unseen planets (Eh] 12004]) . 
However, such models tend to produce discrete rings, 
as around HR 4796A and HD 141569A, rather than 
the radially dimming SBP of the AU Mic disk. Similar 
changes in the power-law index have also been seen in 
the SBPs of other resolved circumstellar d isks: Pic 
llArtvmowicz. Paresce. fc Burro w's! 119901 iHeap et al. 
20001) . TW Hva l)Krist et all 120001: iWeinbereer et al. 
2002) and HD 100546 l|Pantin. Waelkens. fc Lagagf 
20001 lAugereau et al.ll200ll) . A different mechanism, not 
necessarily involving planets, may be at play in these 
systems. 

The TW Hya and HD 100546 circumstellar disks are 
gas-rich and have large mid-plane optic al depths, and 
hence are very much unlike the gas-poor (Ro berge et al.l 
2004), optically thin AU Mic debris disk. However, a 
comparison with Pic is particularly illuminating, be- 
cause of the similar viewing geometry of the two sys- 
tems, and their identical ages. In the remainder of this 
section we seek a common disk architecture that can self- 
consistently account for the broken power-law morphol- 
ogy of the SBPs of these two debris disks. 

The SBP of Pic exhibits a break at 
5-6" (100-120 AU) from the s tar (e.g. 
i Golimowski. Durrance. fc Clampinl 119931 IHeap et all 
2000) , with similar values (from -1 to -4) of the power- 
law index on either side of the break as in the SBP 
of AU Mic. From i-Tz-band (0.21/Ltm) A O observations 
resolv ing the Pic disk over 1.5"-6", iMouillet et al.l 
( 1997) observe the break at a somewhat smaller radius, 
4-4.5" (75-85 AU), with a smaller change in the SBP 
power-law index (from -1 to -3). This may be the inner 
edge of the SBP break observed in the visible, or may 
indicate a wavelength-dependence of the Pic break 
radius similar to the one potentially seen in the AU Mic 
disk (Section l5.1|) . However, because the IMouillet et all 
iO-band data do not extend beyond the optical break 



radius (6"), and because the change in the power-law 
index observed at Ki does not span the full range of 
power-law indices inferred from optical imaging, this 
d ata set will not be c onside red further. 

lArtvmowicz et all l)198 9) model the break in optical 
SBP of the Pic disk using two different power laws 
for the number density of dust particles in the disk, 
for radii less than or g reater than 100 AU, respectively. 
iBackman et alJ l)1992j) consider the possibility of differ- 
ing grain sizes, in addition. While either model may cor- 
rectly describe the architecture of the Pic disk, both 
are purely phenomenological, as they do not model the 
physics behind the discontinuity in the disk. Based on 
the apparent homology between these two debris disks, 
we believe that a plausible two-component model should 
be able to explain both systems self-consistently. In light 
of this, several physical scenarios from the subsequent 
literature are considered below. We find that none of 
them offer a unique explanation, and propose separate 
hypotheses in Sections 15 . 2 . 31 and 15 . 2 . 41 

5.2.1. Ice or Comet Evaporation 

dB^gkm^^ (I1992D and 
mtin. Lagage. fc Artvmowiczl l)1997[) suggest that 
the discontinuity in the SBP of Pic may correspond 
to the location of the "ice boundary" in the disk: all 
dust particles at separations >100 AU are covered with 
ice, while at shorter separations some may not be. This 
results in a deficiency of highly reflective particles in the 
inner regions, creating a shallower power-law index for 
the scattered light profile. In an optically thin disk the 
radius of the ice boundary should scale as the square root 
of the stellar luminosity. Adopting £aumic = 0.13L© 
(Tabled, and Lpp- lc = 8.7L Q , and assuming identical 
grain chemistry in the two systems, we find that the 
corresponding boundary around AU Mic should scale 
down to a radius of 13-15 AU; too close to account for 
the break at 33 AU. 

An alternative hypothesis, involving dust extrac- 
tion through evaporation of gas from a reser- 
voir of cometary bodies around Pic is pro- 
posed by iLecavelier des Etangs. Vidal-Madiar. fc Ferletl 
(1996|. However, the radial distance of this evaporating 
reservoir should scale in the same manner as that of the 
ice boundary. Hence, neither of these two hypotheses can 
be applied simultaneously to AU Mic and Pic. 

5.2.2. A Belt of Parent Bodies 

A reservoir of parent bodies at a discrete range of sep- 
arations from AU Mic could explain the kink in the SBP. 
iGor'kavvi et alJ l)1997[) calculate that the Main asteroid 
belt in the Solar System should produce a break in the 
power-law index of the number density distribution of in- 
terplanetary grains from —1.3 to —6.4 at 0.5-3.0 AU from 
the Sun. These predictions are consistent with empiri- 
cal data from radar meteor s and from impact detectors 
on spacecraft i|Divinei n"993). The inner edge of this belt 
of asteroids in the AU Mic and Pic systems would be 
at the location of the kinks, at ~33 AU and ~110 AU 
from the stars, respectively. By continuing the analogy 
with the Solar System, such belts of parent bodies would 
likely need to be maintained in a discrete range of or- 
bits through mean motion resonances with planets (e.g., 
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iLiou fe ZooklllQQ^ . and references therein). That is, this 
particular scenario may provide further indication for the 
existence of planets in the two disks, in addition to the 
evidence arising from the their clumpy structure. How- 
ever, this model is poorly constrained, as we are free to 
invoke a belt of parent bodies at any distance from either 
star. 

5.2.3. Collisional Evolution 

For two stars of the same age, the disk around the 
more massive star is expected to be colhsionally evolved 
out to a greater radius, because of the inverse scaling of 
the orbital period (and hence, collision frequency) with 
orbital radius and stellar mass. The collisional time scale 
for particles of mean size a on a circular orbit of radius 
r from a star of mass M* is 

p r l/2 r 7+l/2 

<coU 47r 2 a 2 rn(r) * M* /2 7i r-T M* 1/2 n ' 

where P is the orbital period, and no is the normalization 
constant for the number density distribution, that we 
presume scales as Md us t- Given the approximate ratios 
of the stellar masses of AU Mic and 8 Pic (0.28, assuming 
M/3 Pi c = 1.8M Q ), and of their circumstellar dust masses, 
(0.2 ± 0.1, where the mass of the 8 Pic disk was taken as 
the av erag e of the estimates f rom Shc ret. Dent, fc WvattJ 
( 2004) and lDent etUl l|200(y i and assuming 33 AU and 
110 AU as the radii of the kinks in the corresponding 
SBPs, we find that for 7 = 1.5^54 the collisional time 
scales at the respective separations around the two stars 
are equal. This would imply SBPs decreasing approx- 
imately as r _7_1 = r~ 2,5 , which is within the ranges 
found in SectionOfor AU Mic and in lHeap et all ( 2000) 
for 8 Pic, and agrees with the fit to the mean AU Mic 
ff-band SBP (power-law index of -2.3 ± 0.2; Figure gj. 
Therefore, in this scenario the two disk systems scale cor- 
rectly within the errors, indicating that whatever process 
we are observing may scale with the mean time between 
inter-particle collisions. Although the collisional time- 
scale at the location of the break in the AU Mic SBP 
is considerably shorter (1-3 Myr) than the age of the 
star (8-20 Myr), other, slower processes in the disk, e.g., 
grain growth, may scale linearly with the time between 
particle collisions. As noted in Section 15.11 grain growth 
could also explain the observed dependence of grain size 
with orbital radius. 

5.2.4. Poynting- Robertson Drag 

The observed change in the power-law index of the 
SBP may be a reflection of the finite lifetimes of sub- 
micron grains in the inner disk. Having ruled out ra- 
diation pressure as a dominant force on grains larger 
than 0.14/im (Section 15.11) . we propose a hypothe- 
sis based on drag forces for grain removal, in par- 
ticular P-R drag. Although corpuscular drag may 
dominate the dynamics of dust arou nd M stars (e.g., 
iFleming. Schmitt. fc Giampapal 119951) , the strength of 
stellar winds from M dwarfs remains largely unknown. 
We therefore ignore corpuscular drag in the follow ing 
analysis (though see iPlavchan. Jura, fc Lipscvll2005l for 
a discussion of the role of corpuscular drag in the AU Mic 
disk), and consider only P-R drag. As long as the mag- 
nitude of the corpuscular drag force around AU Mic is 



not much greater than the magnitude of the P-R drag 
force, the conclusions remain unchanged. 

If P-R drag was responsible for the depletion of micron- 
sized grains in the inner disks of AU Mic and 8 Pic, then 
the P-R lifetime (ippj of the smallest grains a m j n should 
be constant as a function of disk radius r, and should 
equal the age of the stars (t age )- This can be inferred 
from the expres sion for the P-R lifetime of a particle of 
size a (see, e.g..lBurns. Lamv. fc Soterlll979|k 

where, p is the mean grain density (2.5 g cm -3 for sili- 
cates), and c is the speed of light. Based on the assump- 
tion that P-R drag is the dominant removal mechanism 
for grains larger than the blow-out size (0.14/im), the P- 
R lifetime of the smallest grains is ipR(a m i n ,7') = iage = 
const. Note, however, that the size of the smallest grains, 
Omin, is not a constant, but varies as a mm cx r~ 2 . 

For AU Mic we found a m - m > 0.3/im at 50-60 AU, and 
a m in > 1.6/xm at 17-20 AU ( Section l5"T|) . We obtain 
i PR (0.3^m, 50AU) = 9.9 Myr and t PR .(1.6/mi, 20AU) = 
8.4 Myr. We do not have information about the change 
m Gniin ^S a function of radius i n 8 Pic. We only no t e that 
from mid-IR and visual images . lArtvmowicz et all l)1989fl 
find that they require "few-micron-sized" silicate grains 2 
to model the scattered light at >6.0" (> 1 15 AU) from 
the star. Assuming a minimum grain size of 3/im, we 
find t PR (3pm, 115AU) = 7.8 Myr around 8 Pic. Given 
the uncertainty in a m i n , the obtained P-R time-scales are 
not constrained to better than a factor of 1.5-2. Never- 
theless, they are remarkably similar and agree well with 
the ages of AU Mic and 8 Pic. 

5.2.5. Summary of Proposed Scenarios 

We find that dynamical scenarios based on collisional 
evolution or P-R drag in debris disks offer simpler and 
more self-consistent accounts of the homology between 
the SBPs of the AU Mic and 8 Pic debris disks, com- 
pared to scenarios relying on ice/grain evaporation or 
belts of orbiting parent bodies. Moreover, both our hy- 
potheses can account for the inferred decrease in the min- 
imum grain size with increasing separation from AU Mic. 
We therefore conclude that both are plausible. Given 
the similarities in their predictions, we do not single out 
which one of them is more likely, but defer that analysis 
to a more detailed theoretical work. Regardless of which 
of the two processes is found to be dominant, we can 
confidently claim that optically thin circumstellar disks 
exhibiting breaks in their SBPs are observed in transition 
between a primordial and a debris state. 

6. CONCLUSION 

We have used AO ff-band observations of scattered 
light to probe the morphology of the debris disk at 17- 
60 AU from the young nearby M dwarf AU Mic. We 
find that the disk is within ~1° of edge-on, and that 
it exhibits a number of morphological peculiarities: ra- 
dial asymmetry, spatially resolved clumps, and a change 
in the power-law index of the surface brightness profile 

2 A second solution involving 1-15/^m ice grains is found to be 
equally plausible. However, its grain properties differ widely from 
the ones adopted for the AU Mic circumstellar dust in this paper. 
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near 33 A U. The obse rved morphology agrees with that 
reported in lLiiJ (|20f)4T > . and is suggestive of the existence 
of planetary perturbers in the disk. No planets are de- 
tected down to lM Jup at >20 AU. 

We use a Monte Carlo three-dimensional dust disk 
model to constrain the overall disk parameters, by op- 
timizing them against the AU Mic SED, and near-IR 
and optical scattered light images of the disk (this pa- 
pen lKalas et alJl2004|) . The combined use of SED data 
tracing the thermal emission from large grains, and of 
imaging data tracing grain properties, allows us to break 
several important degeneracies in models of circumstel- 
lar disks that cannot be resolved using only one of the 
two data sets. From the SE D we confirm the previ- 
ously inferred (jLiu et alJl2004|) circumstellar dust mass 
of 0.011M®, and from the properties of the scattered 
light we infer that the debris particles are > 0.5i ; 2 jUm 
in size at 50-60 AU from AU Mic. We find tentative 
evidence for a maximum grain size of 300±2MA«n. How- 
ever, since the data are not sensitive to particles 3>1 mm, 
the result is consistent with no upper limit on the grain 
size. Assuming a single dust size and density distribu- 
tion, we estimate that the radius of the inner disk clearing 
is 10 AU. However, smaller (~1 AU) gap sizes cannot be 
ruled out if a shallower density profile (as observed in- 
wards of 33 AU) and/or larger grains in the inner disk (as 
evidenced from the optical-near-IR color of the disk) are 
adopted. We attribute the lack of sub-micron particles in 
the inner disk either to grain growth, or destruction by 
P-R and/or corpuscular drag (for grains >0.14/im), or to 
blow-out by radiation pressure (for grains <0.14/im). All 
of these mechanisms can explain the increase in relative 
density of small grains with increasing radius in the disk. 

The MC3D model can account for the overall disk pro- 
file and colors to first order, although our one-component 
model fails to reproduce higher-order effects, such as the 
change in power-law index of the SBP. We have discussed 
a number of scenarios which may be capable of reproduc- 
ing such a change in the joint context of the AU Mic and 
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P Pic debris disks. We find that models dividing the disk 
into two separate components, with different grain distri- 
bution and/or composition, represent the combined SED 
and imaging data best. In particular, classes of models 
which scale with the collision or the P-R time-scale are 
most likely to explain both debris disks self-consistently. 

Future high-dynamic range imaging observations prob- 
ing closer to AU Mic (e.g., with nulling interferometry in 
the mid-IR) will further narrow down the architecture 
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vations with Spitzer could better constrain the SED of 
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tures. These could trace small amounts of dust and gas 
in the inner disk, even if no continuum excess is seen at 
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AO Imaging of AU Mic 




Fig. 1. — //-band images of AU Mic. (a) A median-combined image of three sky-subtracted 54.3 sec exposures with the 0.6" coronagraph. 
The star is visible through the semi-transparent coronagraph. The disk is discernible as a faint pair of diametrically opposite spikes along 
the SE-NW direction (traced by the two pairs of parallel lines), (b) An image of the PSF star, HD 195720, with the same coro nagr aph, 
scaled to the intensity of the AU Mic image in panel (a). An arrow points to the location of a faint projected companion (Section l3.3l , (c) 
The final median-combined image of AU Mic obtained from PSF-subtracted images with 0.6", 1.0" and 2.0" spot diameters. All surface 
brightness photometry is performed on this image, (d). Same as in (c), but with an overlaid digital mask covering the star and the telescope 
diffraction pattern. The residual noise in the central region (3.4"-diameter circle) is greater than the surface brightness of the edge-on disk. 
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FIG. 2. — Median ratios /(r) of the radial profiles of AU Mic and HD 195720 (the PSF star) at H band for the images taken with the 
various coronagraphs. The vertical dashed line shows the edge of the largest coronagraph. Our PSF images were multiplied by f(r) before 
subtracting them from the corresponding AU Mic images. For any given coronagraphic spot size, f(r) varies by less than 15% in the range 
r = 1 - 8". 
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Fig. 3. — 12.0" X 4.0" images of the AU Mic disk, with the SE arm oriented horizontally, (a) Locations of the photometry regions for 
measuring the disk surface brightness superimposed on the final, masked, image of AU Mic. The regions are 0.16" X 0.48", and sample 
1-2 FWHMs of the disk thickness. The circular mask is 2" in radius. The crossed arrows mark the location of AU Mic. The contour levels 
trace surface brightness from 17.5 mag arcsec -2 to 14.3 mag arcsec - 2 in s teps o f 0.8 mag arcsec -2 . (b) Small-scale structure in the AU Mic 
disk. The capital letters correspond to sub-structures identified by ILiul 120041) . The bar at the bottom is 10" (100 AU) long, extending 
from —50 AU to +50 AU along the disk plane, with tick marks every 10 AU. To enhance the appearance of the clumps in the disk, we 
have multiplied the pixel values by the square of the distance from the star. The contour levels follow a squared intensity scale, (c) The 
preferred scattered light model of the AU Mic disk at H-band, created using the MC3D code (Section 4;Table[I}. The same software mask 
as in the other two panels has been applied. The contour levels follow the same spacing as in panel (a). No background noise is added, 
though Poisson-noise "clumps" due to low signal-to-noise of the model can be seen. These do not represent discrete physical structures. 



14 



Metchev et al. 



14 



O 
CD 
in 
O 
u 
cd 



!00 

Cti 



w 
w 
cd 
c 

• r— I 

QJ 
O 
Cti 

a 

m 



16 



18 



20 



22 





H band, NW arm 
H band, SE arm 
J? band (Kalas et al. 2004) 
fits to 17-33 AU and 33-60 AU data 
preferred model fit over 17 — 60 AU 



20 



30 40 
Separation (AU) 



50 60 70 



Fig. 4. — H-band surface brightness profiles of the NW (upward pointing solid triangles) and the SE (downward pointing open triangles) 
arms of the AU Mic disk. A gradual flattening of the SBPs of both arms is observed inwards of 30—40 AU. The dashed lines represent the 
power-law fits to t he mean SBP at 17 —33 AU (index of —1.2 ±0.2), and 33-60 AU (index of —4.0 ±0.6). The solid line represents the mean 
ij-band SBP from Kalas et al. (2004) with a power-law index of —3.75. The dot-dashed line shows are preferred model with a power-law 
index -2.2, matching that (-2.3 ± 0.2) of the mean SBP over 17-60 AU. 



AO Imaging of AU Mic 



15 



6 



<X5 



8 



10 



12 



14 



16 



~\ i r 



"i i r 



3 M Jup at 10 Mw 



3 M Jup at 30 Myr 



1 M Jup at 10 Myr 



1 M Jup at 30 Myr 




J I L 



J I L 



10 



50 



100 



Distance (AU) 



Fig. 5. — H-band 5<r detection limits for companions to AU Mic. The solid line delineates the limits in regions away from the disk. 
The break at 10 AU (1") corresponds to the edge of the largest (2"-diameter) coronagraphic spot. To determine the detection limits at 
separations <1", we used only the series of images taken with the 0.6" spot, constituting a third of the total exposure time. Thus, our 
sensitivity at small separations is somewhat worse than what the extrapolation from distances >10 AU would predict. The triangle symbols 
trace the poorer sensitivity to point sources in the plane of the disk. Limited experimentation with planting artificial sourc es in the image 
confir med these detection limits. The dashed lines indicate the expected contrast for 1 and 3Afj up planets around AU Mic (Burrows ct al. 
1997) for system ages of 10 and 30 Myr. 
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Fig. 6. — Examples of degeneracies in the SED. Larger inner gap radii (a) are degenerate with larger minimum grain sizes (b) and shallower 
number density distributions (c). Various combinations of these parameters can produce roughly the same SED. The normalizations of the 
emitted mid-IR flux from the disk are different among the three panels. 
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Fig. 7. — SED of AU Mic (data points from lLiu et al.1120011. and references therein). The photophere is fit by a 3600 K NextGen model 
(Hauschildt ct al. 1999), and the circumstellar excess emission is fit using the MC3D code (Section 4) with model parameters listed in 
Tabled 
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TABLE 1 

Preferred model parameters for the AU Mic system. 



Parameter 


Value 


Range 


Stellar luminosity L„ 


O.13Z/0 


fixed 


Stellar radius 


0.93-Rq 


fixed 


Stellar temperature T c g 


3600 K 


fixed 


Dust mass Md us t 


o.oiiMe 


0.008-0.016Mg) 


Dust size distribution 


-3.5 


fixed 


Number density profile 


r -i.o 


-0.2 to -3.0, fixed 1 " 


Outer disk radius r ou t 


1000 AU 


fixed 


Inner disk radius n n 


10 AU 


1-10 AU 


Inclination angle i 


89° 


> 89° 


Scale height H 


0.8 AU at r = 40 AU 


<1.0 AU at r = 40 AU 


Opening angle 2<5 


0° 


< 6° 


Minimum grain size a m ; n 


0.5^m 


0.3-1/im 


Maximum grain size a max 


300Atm 


100-lOOOAtmt 



1 Treated as a fixed parameter during the MC3D model fitting. The listed 
range corresponds to the range of fits to the SBP over 17-60 AU. 

*The SED data are not sensitive to emission from grains >1000^tm in size, 
so we have not run models with a max > lOOO/um. 



